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ABSTRACT 



We use a simple two-layer outer gap model, whose accelerator consists of 
a primary region and a screening region, to discuss 7-ray spectrum of mature 
pulsars detected by Fermi. By solving the Poisson equation with an assumed 
simple step-function distributions of the charge density in these two regions, the 
distribution of the electric field and the curvature radiation process of the ac- 
celerated particles can be calculated. In the our model, the properties of the 
phase-averaged spectrum can be completely specified by three gap parameters, 
i.e. the fractional gap size in the outer magnetosphere, the gap current in the pri- 
mary region and the gap size ratio between the primary region and the total gap 
size. We discuss how these parameters affect the spectral properties. We argue 
that although the radiation mechanism in the outer gap is curvature radiation 
process, the observed gamma-ray spectrum can substantially deviate from the 
simple curvature spectrum because the overall spectrum consists of two compo- 
nents, i.e. the primary region and screening region. In some pulsars the radiation 
from the screening region is so strong that the photon index from lOOMeV to sev- 
eral GeV can be as fiat as ~ 2. We show the fitting fractional gap thickness of 
the canonical pulsars increases with the spin down age. We find that the total 
gap current is about 50 % of the Goldreich- Julian value and the thickness of the 
screening region is a few percent of the total gap thickness. We also find that the 
predicted 7-ray luminosity is less dependent on the spin down power {Lgd) for 
the pulsars with L^^ > 10'^^ erg/s, while the 7-ray luminosity decreases with the 
spin down power for the pulsars with Lsd ^ 10^^ erg/s. This relation may imply 
that the major gap closure mechanism is photon-photon pair-creation process for 
the pulsars with Lg^ ^ 10^^ erg/s, while the magnetic pair-creation process for 
the pulsars with Lgd ^ 10'^^ erg/s. 
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1. Introduction 

The particle acceleration and high-energy 7-ray radiation process in the pulsar 
magnetosphere has been studied with polar cap model (Ruderman & Sutherland 1975; 
Daugherty & Harding 1982, 1996), the slot gap model (Arons 1983; Muslimov & Harding 
2003; Harding et al. 2008) and the outer gap model (Cheng, Ho & Ruderman 1986a,b, 
hereafter CHRa,b; Hirotani 2008; Takata, Wang & Cheng 2010). All models have assumed 
that the charged particles are accelerated by the electric field along the magnetic field lines, 
and that the acceleration region arises in the charge deficit region from the Goldreich- Julian 
charge density (Goldreich & Julian 1969). The polar cap model assumes the acceleration 
region near the stellar surface, while the outer gap model assumes a strong acceleration 
in outer magnetosphere. The slot gap model proposes that the accelerating electric field 
near the rim of the polar cap can not be completely screened out and the particles are 
continuously accelerated up to high altitude along the magnetic field lines. 

The recent 7-ray instruments have started to break the bottleneck of study of the 
high-energy radiation from the pulsar. In particular, the Fermi LAT has detected 7-ray 
emissions from 46 pulsars just in one year observation, including 21 radio-selected pulsars, 
16 7-ray selected pulsars, the Geminga and 8 millisecond pulsars (Abdo et al. 2010a, 
2009a,b), while the pulsed radio emissions are reported from several 7-ray selected pulsars 
(Camilo et. al. 2009). Romani & Watters (2010) study the pulse profiles of canonical 
pulsars observed by Fermi LAT. They compute the pulse profiles predicted by the outer 
gap model and slot gap model with the magnetic field, which is composed of the vacuum 
dipole field plus the current-induced field. They argued that the outer gap geometry is 
statistically more consistent with the observations than the slot gap geometry. On the other 
hand, the present population study has argued that both outer gap and slot gap model 
may not explain all features of Fermi pulsars; for example, the both model predicts too few 
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young pulsars compared with the Fermi pulsars (c.f. Grenier & Harding 2010). Venter, 
Harding & Guillemot (2009) fit the pulse profiles of the Fermi detected millisecond pulsars 
with the geometries predicted by the different emission models. However, they found that 
the pulse profiles of two out of eight millisecond pulsars cannot be fitted by either the 
geometries with the outer gap or the caustic models. They proposed a pair-starved polar 
cap model, in which the multiplicity of the pairs is not high enough to completely screen 
the electric field above the polar cap, and the particles are continuously accelerated up to 
high altitude over the entire open field line region. 

The observed spectra also allow us to study the site of the 7-ray emissions in the 
pulsar magnetosphere. In the first Fermi pulsar catalog (Abdo et al. 2010a), the spectral 
fits have been done assuming exponential cut-offs. For example, the Fermi LAT measured 
a detail spectral properties of the Vela pulsar (Abdo et al. 2009c, 2010b). It is found that 
the observed spectrum is well fitted by a power law plus hyper-exponential cutoff spectral 
model of the form dN/dE oc E~"- e'XY)[—{E / Ec)'^] with h <1 and Ec ~ 1.5 GeV. This implies 
that an emission process in the outer magnetosphere is more favored than that near the 
stellar surface, which predicts a super exponential cut-off (6 > 1) due to the magnetic 
pair-creation prosses of the 7-rays. Furthermore, the detection of the radiation above 
25 GeV bands associated with the Crab pulsar by the MAGIC telescope has also implied 
the emission process in the outer magnetosphere (Aliu et al. 2008). 

There is a common idea at least in different outer gap models. Charged deficit regions, 
in which electric field along magnetic field lines is not zero, exist in the outer magnetospheric 
regions. The size of outer gap cannot occupy the entire open field line region, because 
electrons/positrons in the gap are accelerated to extremely relativistic energies and each of 
these charged particles radiates a very large number of GeV-photons due to the curvature 
of the magnetic field. Only a small fraction of curvature photons converted into pairs 
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can limit the size of the outer gap (e.g. CHRa,b; Zhang & Cheng 1997). Cheng, Ho & 
Ruderman (1986a) have argued that a reahstic dynamic gap should consist of two regions, 
i.e. a primary acceleration region of lower part of the gap, where seed current can produce 
multi-GeV gamma-ray photons, and a screening region of upper part of the gap, where 
the created electron/positron pairs start to screen out the gap electric field (cf. Fig. 12c of 
CHRa). More detail numerical studies (e.g. Hirotani 2006; Takata et al. 2006) expect that 
a screening region must exist and in fact most gap current is located in this region. 

In this paper we assume a simple two-layer gap structure, which consists of a primary 
acceleration region and a screening region. The electric field along the magnetic field lines 
in both regions are not zero but the exact distribution depends on the charge distribution in 
these two regions. Charged particles accelerated in these regions radiate different curvature 
spectrum due to the position dependent electric field. The combined spectrum depends 
sensitively on the charge distribution in these two acceleration regions. The first purpose 
of this paper is that we explain the photon indies of phase-averaged spectrum of mature 
pulsars by applying our two-layer outer gap model. Observationally the phase-averaged 
spectrum above 100 MeV bands of Fermi 7-ray pulsars can be fit by a single power law 
plus exponential cut-off form (cf. Abdo et al. 2010a). It is found that the photon index 
of the power law component distributes between p ~ 0.5 — 2.5 and the distribution has a 
peak around p ~ 1.3 (Abdo et al. 2010a). It is virtually impossible to explain the very 
soft spectrum with a photon index p > 2 with a simple curvature emission model by 
adjusting the curvature radius of a single layer. We will argue that the two-layer model 
can predict the photon index as flat as p ~ 2 due to the overall spectrum consists of the 
two components, i.e. the curvature radiations from screening region and from the primary 
region. 

The second purpose of this paper is to study the relation between the 7-ray luminosity 
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and the spin down luminosity. The increases of the 7-ray pulsars allow us to perform a 
more detail statistical study. For example, the results of the Ferm will modify the relation 
between the 7-ray luminosity (L^) and the spin down power (Ls^), for which L^ oc L^^ 
with /3 ~ 0.5 were predicted by EGRET 7-ray pulsars (Thompson 2004). However, the 
luminosity L^ calculated from the observed flux include ambiguity due to the uncertainties 
of the distance to the pulsars and of the beaming factor. In this paper, therefore, we 
theoretically predict 7-ray luminosity by fitting the properties of the phase-averaged 
spectrum. This will be important to know the true relation between the 7-ray luminosity 
and the spin down power. 

This paper is organized as follows. In section 2, we present a simple two dimensional 
gap model, which ignores the gap structure in the azimuthal direction. By assuming 
simple step function of the charge distributions in these two acceleration region, we can 
solve the Poisson equation analytically and obtain the electric field everywhere inside 
the gap. We calculate the curvature radiation to compare the observed phase-averaged 
spectrum. In section 3 we discuss how the model parameters affect the spectral properties. 
We also show the fitting results for the phase-averaged spectrum of the all Fermi 
detected gamma-ray pulsars except the Crab-like pulsars, whose radiation spectra could be 
synchrotron-self- Compton process instead of curvature radiation process. In section 4 we 
discuss the implications of the fitted parameters and deduced parameters. We present a 
brief conclusion in section 5. 
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2. Theoretical Model 

2.1. Gap structure with two-layer 

In addition to pulsar parameters, i.e. rotation period and dipole magnetic field, there 
are many factors to affect the observed detailed spectrum of individual pulsar, e.g. the 
phase-resolved spectrum, which depends on the magnetic field structure, the gap structure, 
the inclination angle, viewing angle etc. However the phase-averaged spectrum is less 
sensitive to these factors. In fact it can tell us the most crucial factors, which govern the 
general properties in phase-averaged spectra of all pulsars. It is the aim of this paper to 
use a simple gap model to describe the acceleration and radiation in the outer gap, and 
try to identify the key parameters to affect the spectral properties of pulsars. We consider 
the outer gap accelerator in the magnetic meridian, which includes the rotation axis and 
the magnetic axis. We divide the outer gap accelerator into two-layers, namely the main 
acceleration region and the screening region, in the trans-field direction in the magnetic 
meridian. Fig. [1] illustrates our gap structure with two-layer. 

We express the electric potential in the gap as $ = $' + $0) where $o is the co- rotating 
potential, which satisfies V^^o = —'^'^Pgj with pcj being the Goldreich-Julain charge 
density. In addition, the potential $' is the so called non co-rotating potential, which 
represents the deviation from co-rotating potential and generates the accelerating electric 
field. Using the Poisson equation V^*^ = ~47rp, and assuming the derivative of the potential 
field in the azimuthal direction is much smaller than those in the poloidal plane, we express 
the Poisson equation of $' in the simple 2-dimensional geometry as 

j^ + |.)*'--Mp-.o,,). (1) 

where the coordinates x and z are distance along and perpendicular to the magnetic field 
lines, respectively. Moreover, we assume that (1) the derivative of the potential field in the 



trans-field (z) direction is larger than that along (x-direction) the magnetic field line and (2) 
the variation of the Goldreich- Julian charge density in the trans-field direction is negligible 
compared with that along the magnetic field line. Although these approximations will be 
crude for the pulsars that have the gap size in the trans-field direction comparable to the 
size of the magnetosphere, we may apply these approximations for all pulsars because the 
typical strength of the electric field in the gap might be described by the solution of the first 
order approximation. By ignoring the derivative of the potential field along the magnetic 
field line, we reduce the two-dimensional Poisson equation ([1]) to one-dimensional form that 

— $'(x, z) = -47r[p(x, z) - Pgj{x)], (2) 

where we approximately express the Goldreich- Julain charge density as Pgj{x) ~ ~j^ 
(cf. Eqs. 3.2 and 3.3 of CHRa), where Q is the angular velocity, B is the magnetic field 
strength, s is the curvature radius of the magnetic field lines and s ~ Ric with Ric = c/Q 
being the light radius. 

To solve the equation (|5]) we impose the boundary conditions on the lower and upper 
boundaries. In this study, we consider that the lower boundary is located on the last-open 
field line, and the upper boundary is located on a fixed magnetic field line. For the lower 
boundary, the condition that 

$'(x,^ = 0) = (3) 

is imposed. On the upper boundary, we impose the gap closure conditions that 

$'(a;, z = h2) = and E'^{x, z = hi) = 0, (4) 

where E'^{x, z) = d^' /dz, and h2{x) is the gap thickness measured from the last-open field 
line. The gap closure conditions ensure that the total potential (co-rotational potential -|- 
non co-rotational potential) field in the gap is continuously connected to the co-rotational 
potential field outside the gap. We note that three boundary conditions are imposed for 
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the two order differential equation, implying the location of the upper boundary can not be 
arbitrary chosen, because both Dirichlet-type and Neumann-type conditions are imposed on 
the upper boundary. In fact, the condition E'j_{x, z = /i2) = connects the charge density 
distribution and the gap thickness h2{x) 

As illustrated in Fig. [H we divide the gap into two regions, i.e. the primary acceleration 
region and the screening region. The main acceleration region is expanding between the 
last-open field line and the height z = hi{x) (region I in Fig. [1^), and the screening region 
is expanding between z = hi{x) and z = h2{x) (region II in Fig. [TK). We note that the 
averaged charge density in the outer gap accelerator must be less than the Goldreich- Julian 
value, that is \p\ < \pgj\i in the primary acceleration region to arise a strong electric field 
along the magnetic field line. In the screening region, on the other hand, \p\ > \pgj\ is 
required to satisfy both boundary conditions in equation (jlj). Physically speaking, the 
charge deficit in the primary region is compensated by the charge excess in the screening 
region. 

As demonstrated by the electrodynamic studies (e.g. figure 7 in Hirotani 2006), (1) 
the charge density rapidly increases in the trans-field direction at the height where the 
pair-creation process becomes important, and (2) the increase of the charge density will be 
saturated around upper boundary because of the screening of the electric field due to the 
pairs. In this paper, therefore, we describe the distribution of the charge density in the 
trans- field direction (z-direction) as a step function form (Fig. [T]d); 

J pi(x), if < 2; < hi{x), 
p{x,z) = < (5) 

[ P2{x), if hi{x) < z < h2{x). 

We define the boundary hi{x) with a magnetic line as well as h2{x), implying h2{x)/hi{x) 
is nearly constant along the magnetic field line. 

Using the boundary conditions that ^'^{z = 0) = and $2(2; = ^^2) = and imposing 
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the continuity of the potential field $^ and d^'.^jdz at the height /ii, we obtain the solution 
of the Poisson equation ([2]) as 

{px[x)- pgj[x))z^ ^Cxz, for < z < /ii(x) 

^ {X,Z) = -27T { 

{p2{x)- pgj{x)}{z^ -hl{x)) + Di{z-h2{x)), for hi{x) < z < hiix) 

(6) 



where 



and 



^ . ^ _ (Pl - PGj)hl{hl - 2/12) - (P2 - pGj){hl - /i2) 
Oil^Xj — 

/i2 
(Pl-p2)/i?- ip2- pGj)hl 



Do(x) 



h. 



The gap closure condition E'j_ = on the upper boundary provides the relation among 
the charge densities (pi, P2) and the gap thickness {hi, /i2) that 

[p2(x) - pGj{x)]hl{x) + (pi(x) - P2{x))hi{xf = 0. (7) 

During the fitting process (section l3.2p . we will find that the charge density in the main 
acceleration region is much smaller than that in the screening region, that is, |pi| << |p2|- 
This implies that the typical charge density in the screening region is describes as 
P2(x) ~ Pgj{x)/[1 — (/ii//i2)^]- Namely, both screening conditions $'(z = /i2) = and 
E'j_{z = /i2) = are satisfied when the charge density in the screening region is proportional 
to the Goldreich- Julian charge density. In the realistic situation, therefore, we expect that 
the typical charge density in the screening region is proportional to the Goldreich- Julian 
charge density. For the charge density, pi, in the main acceleration region, we also put the 
charge density proportional to the Goldreich- Julian charge density for simplicity, pi oc pgj, 
because (1) it is expected that the charge density increases along the magnetic field region 
due to the photon-photon pair-creation process and (2) in fact, the distribution of the 
charge density in the main acceleration region is less important to the Potential field in 
the gap (c.f. equation ([6])) as long as the charge density |pi| is much smaller than |p2| and 
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\Pgj\- We will see in section |3^ the observed spectrum can be fitted with pi ~ O.lpcj and 
P2 > Pgj- Because we consider the case that inclination angle is smaller than a < 90° and 
because the outer gap extends beyond the null charge surface, we rewrite the charge density 

as 

P-Pgj ^ \ -9u if Q<z<hi 
PGJ \ g^, if hi<Z<h2' 

where gi > 0, g2 > (cf. Fig. [It). The accelerating electric field along the magnetic field 
line, -E'l = —d^'/dx is described as 

nB \ -giZ^ + C[z, for < z < hi{x) 

E\\{z) S (9) 

g2{z'^ -hl{x)) + D[{z-h2{x)), for hi{x) < z < h2{x) 



cs 



where 



and 



w/x_ gihijhi - 2/^2) + g2{hi - h2f 



D'.ix) 



/i2 

(gi + g2)hl + g2hl 



h2 

where we used d{Bhl)/dx ~ 0, d{z/h2)/dx ~ 0, d{hi/h2)/dx ~ and ds/dx ~ 0. 

A particular solution of equation (|6]) is obtained by specifying any three parameters 
or their combinations out of four, i.e. /ii, /i2, gi and g2- In the next section, we fit the 
observed phase-averaged spectra with the following combination of model parameters, i.e. 
h2{Ric) / Ric, 1 — fi'i and /ii//i2, which correspond to the fractional gap size, the physical gap 
charge density (or current) divided by Goldreich- Julian value in the primary region and the 
gap size ratio between the primary and screening regions. 

2.2. Curvature radiation spectrum 

The charged particles are accelerated by the electric field along the magnetic field lines 
in the gap, and emit 7-rays via the curvature radiation. The Lorentz factor, 7e, can be 
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estimated from the condition that eE'Jz)c = lcur{z), where lcur{z) = 2e'^c'jf{z)/3s'^ is the 
power of the curvature radiation, as 

"3s2 



7e{z) 



2 e-^lW 



1/4 



(10) 



For a single particle, the spectrum of curvature radiation is described by 

where Ecur{z) = {3/2)hc^^{z)/s is the characteristic energy of the radiated curvature 
photon and F{x) = j°° K^i^{y)dy, where K^j^ is the modified Bessel functions of order 5/3. 
The total curvature radiation spectrum from the outer gap is computed from 

F^uriE,) = I n{x,z)FZ''\z)dV, (12) 

where n{x,z) ~ nGj{x){\ + g{z)) is the number density, ncj^x) = VLB{x)/2iice, and 
dV = 6Aj_{x)dx is the volume of the element, where 6A±{x) is the increment area of 
width 6z. The increment are is calculated as 6A±{x) ~ 2'jtRicSz{Ric)B{Ric)/B{x), where 
Sz{Ric) is the width at the light cylinder and we used the magnetic flux conservation that 
B{x)SA^{x) = B{Ri,)5A^{Ri,) ~ 2tiBi,Ri,5z{Ri,). 

The total flux received at Earth is 
where d is the distance of the pulsar and Afi is the solid angle. 



3. Results 

3.1. Properties of the curvature spectra with the gap structures 

Basically we have three independent fitting parameters, i.e the fractional gap size 
defined by / = h2{Ric)/Ric, the number density (current) in the primary region (1 — gi) and 
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the ratio between the thicknesses of primary region and total gap size /ii//i2- Quahtatively 
the fractional gap size mainly determines the total potential drop in the gap and the 
strength of the accelerating field, and hence it affects the intensity and the cut-off energy 
of the curvature radiation. We will see that the number density {1 — gi) in the primary 
region determines the slope of the curvature spectrum below the cut-off energy and the 
ratio {hi/h2) controls the spectral break in lower energy bands around 100 MeV. In below 
we give the quantitative description how these three parameters affect the spectrum. 

Fig. [2] summarizes the dependency of the curvature spectra on the gap thickness. The 
lines present the spectra with the fractional gap thickness of /=0.8 (solid line), 0.6 (dashed 
line), 0.4 (dotted line) and 0.2 (dashed-dotted line). The results are for 1 — gi = 0.10, 
hi/h2 = 0.933 and the parameters of the Geminga pulsar. We can see in Fig. [2] that the 
spectral properties are sensitive to the gap thickness. Specifically, the intensity increases 
with the gap fractional thickness. This is because the total potential drop in the gap and 
therefore magnitude of the accelerating electric field are proportional to the square of the 
gap thickness ($' oc /^), as equations (EI) and (Q imply. Therefore, the predicted flux of 
the curvature radiation from the gap is approximately proportional to E^F oc $' x A^ oc /^, 
where we used that the total number of the particles (or current) in the gap is proportional 
to N (X f. Fig. [2] also shows that the spectrum becomes hard with increase of the fractional 
gap thickness. From the equations iQ and flTOj) . we can see that the typical energy of the 
curvature photons is proportional to E^ oc {3/2)hc'y^/s oc /'^^^, implying the cut-off energy 
in spectrum increases with the fractional gap thickness as Fig. [2] shows. 

Fig. [3] shows how the charge density, p = (1 — gi)pGJ, of the primary region affects the 
curvature spectrum. The lines represent the curvature spectra for 1 — gi = 0.30 (solid line), 
0.10 (dashed line), 0.05 (dotted line) and 0.01 (dashed-dotted line). The results are for 
/ = 0.76, hi/h2 = 0.933 and the parameters of the Geminga pulsar are used. Qualitatively 
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speaking, the number density in the primary region determines the slope of the spectrum 
below the cut-off energy. We can find in Fig. [3] (e.g. dashed-dotted line) that the spectrum 
is divided into two components, that is, higher and lower energy components. The higher 
energy component {E^ > 1 GeV) comes from the primary region, where the charge density 
is smaller than the Goldreich-Julain charge density, while the other one {E^ < 1 GeV) 
originates from the screening region, where the charge density exceeds the Goldreich-Julain 
charge density. Because the accelerating electric field in the primary region is stronger than 
that in the screening region, the Lorentz factor of the particles and the resultant typical 
energy of the curvature photons are larger in the primary region than in the screening 
region. 

The screening condition, (/i2/^i)^ = 1 +91/92 (c-f- Eq|7]), implies a larger number 
density (corresponding to larger (72) in the screening region is required for a smaller number 
density (corresponding to larger gi) in the primary region. In such a case, if the number 
density of the primary region is much smaller than Goldreich- Julian value, the energy flux 
of the curvature radiation from the primary region becomes small compared with that from 
the screening region, as indicated by the dashed-dotted line in Fig [31 With the present 
model, therefore, number density in the primary region, and resultant ratio of the total 
numbers of particles in the primary and the screening regions mainly determines the slope 
of the spectrum below the cut-off energy. Fig. |3] also shows that the cut-off energy in the 
spectra in GeV energy bands decreases as the charge density in the primary region (or 
1 — Qi) increases. This is because the deviation from the Goldreich- Julian charge density 
becomes small as the charge density increases, implying that the accelerating electric field 
and the resultant typical energy of the curvature photons decrease with increasing of the 
charge density. 

FigJH summarizes the dependency of the properties of the curvature spectrum on the 
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ratio of the thicknesses between the primary region and the screening region, /ii//i2- The 
hnes represent the curvature spectra for /ii//i2 = 0.70 (sohd hne), 0.80 (dashed hne), 0.90 
(dotted hne) and 0.96 (dashed-dotted hne) with / = 0.76, 1 — gi = 0.1 and the parameters 
of the Geminga pulsar. We note that the increase of the ratio /ii//i2 corresponds to the 
decrease of the thickness of the screening region (/i2 — hi) relative to the primary region. 
In the present model, the ratio /ii//i2 determines the energy width between the spectral 
cut-off in GeV band and the spectral break in 100 MeV bands, because the curvature 
spectrum of the primary region has an energy peak at several GeV, while that of the 
screening region shows an energy peak at several hundred MeV. In the dashed-dotted line 
in Fig. m {hi/h2 = 0.96), for example, the spectral cut-off appearing at several GeV is 
a result of the emissions from the primary region, while the spectral break seen around 
200 MeV is produced by the emissions from the screening region. As the ratio decreases, 
on the other hand, the position of the spectral break in lower energy shifts to higher energy 
(e.g. ~500 MeV for /ii//?-2 = 0.90), and the energy width between the cut-off energy and 
the spectral break energy decreases. This is because the thickness of the screening region 
(/i2 — hi) relative to that of the primary region (hi) becomes thick as the ratio decreases, 
implying the flux of the emissions from the screening region increases and its spectrum 
becomes hard. For the ratio hi/h2 = 0.70 (solid line in Fig. H]), the emissions from the 
screening region dominates the emissions from the primary region. 

3.2. Fitting Results 

We fit the phase-averaged spectra of the 42 7-ray pulsars measured by the EGRET 
and Fermi telescopes with the present model. For the observed spectra, we use information 
reported in the first Fermi catalogue (Abdo et al. 2010a), in which the observed data 
were fit with a single power law plus exponential cut-off form. In our fitting, we do not 
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include the Crab- like young pulsars (e.g. the Crab pulsar, PSR Jl 124-5916) because the 
radiation mechanism of the Crab-like pulsars are synchrotron-self-Compton process instead 
of curvature radiation process (Cheng, Ruderman & Zhang 2000; Takata & Chang 2007). 
In the Crab-like pulsars their soft photon density is sufficiently high so that most curvature 
photons from the outer gap will be converted into pairs within the light cylinder and the 
observed gamma-rays resulting from the inverse Compton scattering of the synchrotron 
photons of secondary pairs. This radiation process differs from what we have considered in 
this paper and therefore we will not consider them. 

As discussed in section 13. Ij our fitting parameters are the fractional gap thickness /, 
the charge density in the primary region, 1 — gi, and the ratio /ii//i2- In fact, given observed 
phase-averaged spectrum can be uniquely fit with one set of (/, 1 — gi, hi/h2) with a small 
uncertainty. The model fitting was proceeded as following. First, we deduced the typical 
fractional gap thickness, /, from the observed intensity and the cut-off energy, because 
the fractional gap thickness / greatly affects the intensity and the cut-off energy as Fig. [2] 
shows. For the next step, we fit the spectral slope below the cut-off energy with the charge 
density in the primary region, 1 — gi, which mainly controls the slope of the calculated 
spectrum as Fig. [3] shows. Finally, we determined hi/h2, which controls the spectral width 
between the cut-off energy in several GeV and the spectral break energy in lower energy 
bands. 

Fig. [5] presents the fitting results with the observed data for the 6 canonical pulsars. 
The data points are taken from the EGRET observations (Fierro, 1995) for the Geminga, 
PSRs J057-5226,J1709-4229 and J1952+3252, and from the Fermi observations for the 
Vela (Abdo et al. 2009c and 2010b) and PSR J2021+4206 (Trep et al. 2010). The grey 
strips represent the errors of the photon index, cut-off energy and intensity measured by 
the Fermi observations. The solid lines represent the best fit spectra with the fitting 
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parameters listed in each panel and in Table 1. We use the dashed and dashed-dotted lines 
in the panels of the Gemiga and Vela pulsars to preset how the best fit parameters include 
the uncertainties. The dotted lines are results for the fractional gap thickness (/ = 0.7 for 
the Geminga and 0.145 for the Vela) about 10 % difference than the best fitting values, 
while the dashed-dotted lines are results for 1 — gi (0.13 for the Geminga and 0.1 for the 
Vela) and hi/h2 (0.867 for the Geminga and 0.967 for the Vela) about 10-20 % difference 
than the best fitting values. It is obvious that both dashed and dashed-dotted lines can not 
explain the observed data, implying the fitting parameters include uncertainties of about 
10 % for the Geminga and Vela pulsars. 

Figs [6] and [7] compare fitting results with the Fermi observations for the 28 canonical 
pulsars and 8 millisecond pulsars, respectively. The solid line and dashed line are 
corresponding to the model spectra and the observations, respectively. The grey strips 
represent the errors of the observations. In Figs IMZl "we can see that the present model 
can reproduce well the observations, although there is a small discrepancy between the 
model and the observations around 100 MeV energy bands for some pulsars (such as PSRs 
J1907+06 and J2229+6144). 

In Table 1, we summarize the observed pulsar parameters (second- fifth columns) and 
the fitting parameters (sixth-ninth columns). Fig. |5] plots the fitted fractional gap thickness 
as a function of the spin down age r. We find in Fig. [5] that the fitted fractional gap size 
ffit (sixth column) is between 0.1-0.9 and is bigger than the estimated fractional gap size 
of the Crab pulsar (e.g. CHRb estimated fcrab ~ 0.05). We can see a clear trend in Fig. E] 
that the fractional gap thickness for the canonical pulsars increases with the spin down age. 
This is expected because the potential drop of the pulsar decreases as the spin down age 
increases, implying a larger fractional gap thickness is required to accelerate the particles 
to emit several GeV 7-ray photons via the curvature radiation process. For the millisecond 
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pulsars, it is required more number of sample to argue the relation between the fractional 
gap size and the spin down age. 

In seventh column of Table 1 we can see that the current in the primary region is 
roughly ~ 10% of the Goldreich- Julian current. This is important because if the gap current 
in this region is too large then accelerating electric field is not strong enough to accelerate 
the gap electrons/positrons to sufficiently high Lorentz. Consequently multi-GeV photons 
cannot be generated and pair creation process cannot take place. 

The ratio between the primary region and the total gap size hi/h2 is less than but 
very closed to unity. This indicates that the screening region is very thin and is only a 
few percents of the total gap size. This is expected because although the mean free path 
(A) of photon-photon pair-creation is longer than the light cylinder radius in the case of 
mature pulsars, the thickness of the screening region is substantially reduced by the photon 
multiphcity (A'^) and it can be estimated as 

(cf. equation 5.7 and C21 of CHRa). Taking A ~ l^^Ric, s ^ Ric ^ lO^cm and N^ ~ 10^ 
h2{Ric) — hi{Ric) ~ 10^ cm << h2{Ric) ~ 10*^ — lO^cm. Also the magnetic pair-creation near 
the stellar surface could supply the screening pairs within a thin layer (Takata et al. 2010). 

In order to match the observed gamma-ray flux, we need to introduce one more fltting 
parameter, i.e. AQd"^ (ninth column), which is the product of gamma-ray solid angle and 
the distance to the pulsar. Using the observed distance (fourth column) and the fltting 
AQd"^, we estimate the solid angle AQfu (twelfth column). We can see that the averaged 
value of the deduced solid angle AQfu is order of unity, which is usually assumed. On the 
other hand, some 7-ray pulsars are required a solid angle much larger (> 10) or smaller 
(< 0.1) solid angle than the averaged value. This would reflect the effects of the viewing 
geometry. The fltting solid angle much larger (or smaller) than the mean value may imply 
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that the pulsars are observed with a viewing angle which cuts through the emission region 
where the intensity is much smaller (or larger) than the mean intensity. 

For the observed distance to pulsar, we can see that first of all the error of the observed 
distance is usually quite large and secondly there are some pulsars without the measured 
distance. Therefore if the solid angle of gamma-ray pulsars is the same, which is usually 
assumed to be unity, then we can use the theoretical predicted power and the observed flux 
to predict the distance to the pulsar, which is listed in the last column of Table 1. 

4. Discussion 

In Table 1 we deduce the averaged gap current in units of the Goldreich-Julain value 
(tenth column). In addition to show that the curvature radiation emitted by the accelerated 
charged particles can consistently explain the observed spectrum of mature pulsars detected 
by Fermi, it is interesting to note that from Table 1 that the gap current is roughly 
constant. Specifically, the variation of the gap current is less than 20% and the gap current 
is around 50-60% of the full Goldreich- Julian current. Although it is not easily explain why 
the total current in the gap is close to half of Goldreich- Julian value, this result appears 
qualitatively reasonable because if the current is very closed to the Goldreich- Julian current, 
the electric field will be substantially reduced. Consequently the charged particles cannot 
be accelerated to extremely relativistic to emit multi-GeV photons and hence pair creation 
process cannot occur. In fact we speculate that the gap may operate in a dynamical 
format, namely, the gap can begin with an almost charge starvation state, and then charged 
particles are accelerated to extremely relativistic and emit 10^ — 10^ photons per particle 
across the gap. A small fraction of these photons are converted into pairs, which can almost 
immediately quench the gap, and the current in this stage is almost full Goldreich- Julian 
current. Since particles and photons are all moving in speed of lights, the on/off time scales 
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of these two stages should be very similar and which may be order of Ric/c. After averaging 
over time, it may not be surprised to have a time average gap current about 50% of the 
Goldreich Julian value. If there were good time resolution, this speculation predicts that 
there is a time variability in the 7-ray intensity with a time scale of Ric/c. This may be 
an evidence of the pair-creation process in the pulsar magnetosphere, which maintains the 
global electric circuit of the star. 

The predicted 7-ray luminosity {L^'^^), for each pulsar is listed in eleventh column 
of Table 1, where we use L;^** = fjifLgd with L^^ = {2'jT)^B'^Ry6c^P^ being the spin 
down luminosity. We note that uncertainties of the predicted 7-ray luminosity is small 
because the uncertainties of the fitting parameters is small, say about 10 %, as discussed 
in section 13.21 In Fig. [9l the predicted 7-ray luminosity is plotted as a function of the spin 
down power Lsd- We can see a trend in Fig. [9] that the predicted 7-ray luminosity shows 
less dependency on the spin down power for the pulsars with L^^ > 10^^ erg/s, while L^ 
decreases with the spin down power for the pulsars with L^^ < 10^^ erg/s. This change of 
the dependency of the 7-ray luminosity on the spin down power may be caused by switching 
the gap closure process. In fact, the solid and dashed lines in Fig. |9] represent the relation 
predicted by the outer gap model with different gap closure processes as follows. 

Zhang & Cheng (1997) have argued that when the pulsars cool down, the cooling 
X-rays may not be sufficient to convert the curvature photons emitted by the accelerated 
particles into pairs to restrict the growth of the outer gap. They suggest that however 
the X-rays emitted by the heated polar cap due to the return current can provide a 
self-consistent mechanism to restrict the fractional size of the gap as 

fzc = 0.32PT'b;,'/\ (15) 

where P_i is the rotation period in units of 0.1s. This model predicts the 7-ray luminosity 
is related with the spin down luminosity as Ly oc LJ^ B^2 ~ LJ^ , which is represented by 
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the solid line in FigjQ] This less dependency on the spin down luminosity may explain the 
behavior of the predicted 7-ray luminosity of the pulsars with Lg^ > 10^^ erg/s. 

Recently Takata et al. (2010) argue that since the electric field decreases rapidly from 
the null charge surface to the inner boundary, where the gap electric field must vanish, the 
radiation loss cannot be compensated by the acceleration of the local electric field. When 
this occurs the Lorentz factor of the incoming electrons/positrons is determined by the 
equating the radiation loss time scale and the particle crossing time scale. It is interesting 
to note that the characteristic photon energy of curvature radiation is independent of pulsar 
parameters and is given by mgC^/af ~ lOOMeV, where a/ is the fine structure constant. 

They further argue that these lOOMeV photons can become pairs by magnetic pair 
creation process. These secondary pairs can continue to radiate due to synchrotron radiation 
and the characteristic energy of synchrotron photons is of order of several MeV. The photon 
multiplicity is easily over 10^ per each incoming particle. Such cascade process has also 
been considered before (e.g. Cheng & Zhang 1999). For a simple dipolar field structure, 
all these pairs should move inward and they cannot affect the outer gap. However they 
argue that the existence of strong surface local field (e.g. Ruderman 1991, Arons 1993) has 
been widely suggested. In particular if the field lines near the surface, instead of nearly 
perpendicular to the surface, are bending sideward due to the strong local field. The pairs 
created in these local magnetic field lines can have an angle bigger than 90°, which results 
in an outgoing flow of pairs. In fact it only needs a very tiny fraction (1-10) out of 10^ 
photons creating pairs in these fleld lines, which are sufficient to provide screening in the 
outer gap when they migrate to the outer magnetosphere. 

They estimate the fractional gap thickness when this situation occurs as 

fm = 0.25i^P!f , (16) 

where K oc B^^2-^j is the parameter to characterize the local parameters, e.g. -Bm,i2 and 
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Sr are the local magnetic field in units of lO^^G and the local curvature radius in units of 
lO^cni. From this estimate they predict that the gamma-ray luminosity is related to the 
spin down power as L^ oc LJ^ , which is represented by the dashed-line in Fig [Hi This gap 
closure process may explain the relation between the predicted 7-ray luminosity and the 
spin down power for the pulsars with Lgd < 10^^ erg/s, as Fig. [9] indicates. 

5. Conclusion 

In this paper, we applied the two-layer outer gap model to fit the observed phase- 
averaged spectra of the 42 7-ray pulsars detected by the Fermi telescope. Our gap 
structure consists two parts, which are the primary acceleration region and the screening 
region (Fig. [1]). In the primary acceleration region, the charge density is less than the 
Goldreich- Julian charge density (|p| < |pgj|); while in the screening region, the charge 
density exceeds pcj to screen out the accelerating electric field. Assuming a step function 
of distribution of charge density in the direction perpendicular to the magnetic field lines, 
we solve the Poisson equation to obtain the accelerating electric field in the gap. We fit the 
observed phase-averaged spectrum with the three fitting parameters, that is, the fractional 
gap thickness /, the number density of the primary region (1 — gi), and the ratio of the 
thicknesses of the primary and the secondary regions, /ii//i2- We demonstrated that the gap 
thickness affects the position of the cut-off energy (several GeV) appeared in the curvature 
spectrum, as Fig|2] shows. We also showed that the number density of the primary region 
and the resultant ratio of the number densities between the primary and the screening 
region determine the slope of the spectrum below the cut-off energy (Fig. [3]), and that the 
ratio /11//12 mainly affects the position of the spectral break at 100 MeV bands (FigS]), 
which is caused by the emissions from the screening region. With two-layer mode, the 
observed soft spectrum with the photon index of ~ 2 for the some pulsars can be explained 
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by the overall spectrum consists of the two components, i.e. the curvature radiations from 
screening region and from the primary region. Our fitting results show that the fractional 
gap thickness for the canonical pulsars tends to increases with the spin down age, as FigO 
shows. The observations can be fit by the ratio hi/h2 smaller than but close to unity, 
implying the screening region is only a few percent of total gap thickness. The present 
model predicts the gap current is about 50 % of the Goldreich-Julain value. We found 
that the predicted 7-ray luminosity shows less dependency on the spin down power for the 
pulsars with Lgd ^ 10^^ erg/s, while it decreases with the spin down power for pulsars 
with Lgii < 10^^ erg/s (Fig. [H]). We discussed the relation of the 7-ray luminosity and the 
spin down power with the gap closure mechanisms of Zhang& Cheng (1997), which predict 
Ly oc LJ^ and of Takata et al. (2010), which predicts L^ oc LJ^ . 

While the present simple model can be applied to discuss the observed phase-averaged 
spectra, a three-dimensional model should be required to discuss the observed light curves 
and the phase-resolved spectra. The detail properties of the observed light curves (e.g. 
number of the peak and phase of peak) and the phase-resolved spectra will reflect the 
three-dimensional distributions of the number density of particle and electric field in the 
gap. For example, the Fermi telescope revealed the third peak, whose position depends on 
the energy bands, in the light curve of the Vela pulsar (Abdo et al. 2009d, 2010). The third 
peak in GeV energy bands, which has not been expected by the previous emission models, 
will be understood with a detail three-dimensional analysis. In the subsequent papers, 
therefore, we will extend the present two-dimensional analysis into a three-dimensional one, 
and fit the phase resolved spectra and the energy-dependent light curve for the individual 
pulsar. 

We wish to express our thanks to the referee for his/her insightful comments on the 
manuscript. We thank the useful discussions with H.-K. Chang, K. Hirotani, C.Y. Hui, 
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Fig. 1. — Schematic view of the gap structure. (a)The structure of the outer gap with two- 
layer. The regions / and // are the primary region and the screening region respectively. 
(b)The distribution of the charge density in the trans-field [z) direction in the gap. (c)The 
distribution of the g{z) in z-direction in the gap. (d)The distribution of the $'(2;), and the 
accelerating electric field E'Jz) oc $'(2;). 
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Fig. 2. — The effects of / on the shape of the spectrum of the curvature radiation. The 
results are for 1 — gi = 0.1 and hi/h2 = 0.933. The parameters of Geminga are used. 
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Fig. 3. — The effect of the charge density of the main acceleration region, 1 — gi, on the 
shape of the spectrum, where /11//12 = 0.933 and / = 0.76 are used. The parameters of 
Geminga are used. 
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Fig. 4. — The effect of /ii//i2 on the shape of the spectrum with the l—gi = 0.1 and / = 0.76. 
The parameters of Geminga are used. 
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Fig. 5. — The model results of 6 canonical pulsars. The solid lines represent the best fitting 
model spectra with the fitting parameters listed in each panel. The circles are the observed 
data from the Fermi LAT, which are taken from Abdo et al. (2009d) for the Vela pulsar 
and Trep et al. (2010) for J2021+4026. The boxes are the observed data from EGERT 
(Fierro 1995). The dotted lines in the panels of the Geminga and the Vela pulsars are results 
for the fractional gap thickness (/ = 0.7 for the Geminga and 0.145 for the Vela) about 
10 % difference than the best fitting fractional gap thickness, while the dashed-dotted lines 
are results for 1 — gi (0.13 for the Geminga and 0.1 for the Vela) and hi/h2 (0.867 for the 
Geminga and 0.967 for the Vela) about 10-20 % difference than the best fitting parameters. 
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Fig. 6. — The model results of 28 canonical pulsars. The solid lines are model results, the 
dash lines correspond to the photon indexes, cut off energies and photon fluxes from the 
Fermi catalogue and the grey strips represent the errors of the three(Abdo et al. 2009a). 
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Fig. 7. — The model results of the millisecond pulsars. The solid lines are model results, the 
dash lines correspond to the photon indexes, cut off energies and photon fluxes from the 
Fermi catalogue and the grey strips represent the errors of the three(Abdo et al. 2009a). 
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Fig. 8. — The fitting fractional gap tliickness as a function of the spin down age. The circle 
and triangle represent for the canonical and millisecond pulsars, respectively. 
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Fig. 9. — The predicted 7-ray luminosity as a function of the spin down power. The circle 
and triangle represent for the canonical and millisecond pulsars, respectively. The solid and 
dashed lines represent the models by Zhang & Cheng (1997) and Takata et al. (2010), 
respectively. 
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0.08-1.72 



1.37-32 



9-1-1.54 
^-1.09 

1 c: + 0.32 
J^-J_0.29 

1-1-3.28 
1^-0.64 

,-1-0.13 
^-0.12 

0.36ig:?» 

„ 27+0.28 
'^■^'-0.11 

1 4+1-28 <^ 

0.09-0.55 



0.04 



-1-0.04 
-0.02 



fl.l4 
-0.44 

0.05-0.32 

fO.18 
■0.07 



1.05 



0.17" 



2.5 



-2.45 
-1.00 

0.05-0.30 

hi. 22 



0.25" 



■0.15 



2.12 
2.62 
0.85 
4.24 
2.19 
0.35 
0.35 
2.07 
0.29 
1.40 
1.41 
0.85 
1.48 
1.12 
1.20 
2.66 
0.79 
1.58 
1.27 
0.60 
1.10 
0.85 
1.12 
1.13 



Table 1:: continued 



Name 



P (ms) Bi 



Observed Parameters 
d°''= (kpc) F°5'o={10-*phcm-2s 



') ff. 



Fitting Parameters 
1-gi hi/h2 AQd'^ (kpc^) 



Deduced Parameters 
rjgap Li,'' (1033erg/s) AQj^t 



d{An = 1) 



J1836+5925* 


173 


0.51 


<0.8 


65.6 ±1.8 


0.90 


0.10 


0.940 


0.32 


0.564 


8.748 


> 0.5 


0.56 


J1907+06* 


107 


3.09 




40.25 ± 3.8 


0.26 


0.05 


0.920 


3.81 


0.545 


49.92 




1.95 


J1952+3252 


39.5 


0.48 


2.0 ±0.5 


17.6 ±1.9 


0.22 


0.08 


0.927 


6.6 


0.558 


39.82 


1 fi^+l-28 
l-OO_0.59 


2.57 


J1958+2846* 


290 


7.95 




7.65 ±1.6 


0.38 


0.20 


0.967 


8 


0.607 


19.64 




2.83 


J2021+3651 


104 


3.18 


2 1+2-1 
^■-^-1.0 


67.35 ± 4.4 


0.18 


0.07 


0.927 


0.7938 


0.553 


19.71 


o.i8lg:S 


0.89 


J2021+4026* 


265 


3.84 


1.5 ±0.45 


152.6 ±4.9 


0.50 


0.06 


0.927 


0.28125 


0.548 


14.5 


0.125tg:J^ 


0.53 


J2032+4127* 


143 


1.68 


1.6-3.6 


6 ±2.3 


0.58 


0.30 


0.953 


25.2 


0.658 


51.31 


1.94-9.84 


5.02 


J2043+2740 


96.1 


0.35 


1.80 ±0.54 


2.41 ±0.90 


0.33 


0.12 


0.933 


2.961 


0.575 


2.012 


n q+094 


1.71 


J2229+6114 


51.6 


2.00 


0.8-6.5 


32.6 ±2.2 


0.10 


0.07 


0.940 


2 


0.549 


22.5 


0.05-3.125 


1.41 


J2238+59* 


163 


4.04 




6.8 ±1.8 


0.20 


0.15 


0.967 


3.64 


0.582 


7.224 




1.91 



Millisecond pulsars 



Observed Parameters 



Name 



P (ms) Ba 



jobs 



(kpc) 



F°o''^(10"**phcm- 



') ff. 



Fitting Parameters Deduced Parameters O^ 

1-gi hi/h2 Af7d2 (]^pf.2) ^^^p L^" (lO^^erg/s) AQfu d{An = 1) 



J0030+0451 

J0218+4232 

J0437-4715 

J0613-0200 

J0751+1807 

J1614-2230 

J1744-1134 

J2124-3358 



4.9 


2.27 


0.3 ±0.09 


2.3 


4.14 


2.5-4 


5.8 


2.9 


0.1563 ±0.0013 


3.1 


1.76 


o.48l°;l? 


3.5 


1.5 


0.61-1 


3.2 


1.2 


1.27 ±0.39 


4.1 


1.8 


U.JO /_Q 035 


4.9 


2.4 


o 9C-I-0.25 



5.83 ±0.78 

6.2 ±1.7 
3.65 ±0.84 
3.38 ±0.85 
1.35 ±0.66 

2.89 ±1.2 

4.3 ±1.6 
1.95 ±0.49 



0.60 


0.12 


0.947 


0.36 


0.572 


0.20 


0.05 


0.920 


1.17 


0.545 


0.38 


0.05 


0.927 


0.12215 


0.543 


0.46 


0.08 


0.947 


0.8064 


0.553 


0.56 


0.07 


0.947 


1.62 


0.548 


0.55 


0.10 


0.933 


0.80645 


0.566 


0.33 


0.15 


0.953 


0.15294 


0.585 


0.80 


0.15 


0.953 


2 


0.585 



0.648 
1.92 
0.165 
1.265 
1.054 
0.832 
0.1438 
2.048 



4+4.16 
^-1.63 

0.07312-0.1872 

c+O.OS 



3.5 



+2.39 



4.5 
0.5 
1.2 
32-* 



+5.63 
-3.38 

+0.54 
-0.21 

+0.28 



0.6 

1.08 

0.35 

0.898 

1.272 

0.898 

0.391 

1.414 



Note. — The first column is tlie name of tlie pulsar, tlie pulsar -witli "*" means it is a 7-ray selected pulsar, which 
was detected by Fermi blind search. The second to fifth column are observed parameters: periods, surface magnetic 
fields in units of lO^^G, observed distances and the photon flux respectively. The data of these four columns are from 



the pulsar catalogue of the Fermi LAT(Abdo et al. 2010a). gi, hi/h2, ffit{Ric) and AQ(P are fitting parameters. 



rjgap = -^ = - — ^ '^ fa ^ — ^^ ^^^ average gap current in units of Goldreich- Julian current. Ll^* = f%^Lsd is 



PG.J 



the 7-ray luminosity of the model result. AQfu = AQd'^/d'^f^^, the errors are due to the errors of the distances. The 
d{AQ = 1) is the predicted distance in units of kpc, it is obtained from the fitting parameter AQd"^, when AQ = 1. 
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